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ABSTRACT 

The warm-hot intergalactic medium (WHIM) at temperatures 10 5 - 10 7 K is believed to contain 30-50% of 
the baryons in the local universe. However, all current X-ray detections of the WHIM at redshifts z > are of 
low statistical significance (<3er) and/or controversial. In this work, we aim to establish the detection limits 
of current X-ray observatories and explore requirements for next-generation X-ray telescopes for studying the 
WHIM through X-ray absorption lines. We analyze all available grating observations of Mrk 421 and obtain 
spectra with signal-to-noise ratio (S/N) of ~ 90 and 190 per 50 mA spectral bin from Chandra and XMM- 
Newton observations, respectively. Although these spectra are two of the best ever collected with Chandra and 
XMM-Newton, we cannot confirm the two WHIM systems reported by Nicastro et al. in 2005. Our bootstrap 
simulations indicate that spectra with such high S/N cannot constrain the WHIM with O vn column densities 
Novn ~ 10 15 cm" 2 (corresponding to an equivalent widths of 2.5 mA for a Doppler velocity of 50 km s -1 ) 
at >3er significance level. The simulation results also suggest that it would take > 60 Ms for Chandra and 
140 Ms for XMM-Newton to measure the A/ovn at > 4a from a spectrum of a background QSO with flux of 
~ 0.2 mCrab (1 Crab = 2 x 10~ 8 erg s" 1 cm" 2 at 0.5-2 keV). Future X-ray spectrographs need to be equipped 
with spectral resolution R ~ 4000 and effective area A > 100 cm 2 to accomplish the similar constraints with 
an exposure time of ~ 2 Ms and would require ~ 1 1 Ms to survey the 15 QSOs with flux > 0.2 mCrab along 
which clear intergalactic O vi absorbers have been detected. 



Subject headings: intergalactic medium — X-rays: diffuse background 
Lacertae objects: individual (Markarian 421) 



quasars: absorption lines — BL 



1. INTRODUCTION 

Identifying the "missing baryons" is one of the major tasks 
of modern co smology. Observation s of the microwave back- 
ground (e.g., iKomatsu et al.l l201lh and the big bang nucle- 
osy thesis model combin ed with the measu r ement of deu- 
terium abundance (e.g., iBurles et aD 120011 : lO'Meara et alJ 
2006; iPettini et all [2008) are converging on the cosmologi- 
cal baryon density of fij = 0.0455 ± 0.0028. At high red- 
shift (z > 3) universe, baryons exist primarily in the form of 
cool, photoionized intergalactic medium (IGM) th at is traced 
by Lya absorption forest lines (Rauc het al.l 119971) . In the 
present-day universe, matter detected in forms of photoion- 
ized IGM, stars, galaxies, intracl uster medium, etc ., adds up 
to only ~ 50% of the baryons (IShull et all 1201 lb. leaving 
much of the inventory yet to be found (Fukugita et al.|[l998T; 
Fukug ita & PeebleTI 2004). Cosmological numerical simula- 
tions for large-structure formation indicate that the missing 
baryons are still in the IGM, but at current epoch they have 
been heated by gravitational shocks and galactic feedback to 
temperatures of 1 5 — 1 7 K when they fall into the gravita- 
tional potential wells of the dark ma t ter cosmic web fil aments 
(e.g., ICen & Ostrikerl 119991 120061: iDave et al.l 120011) . The 
IGM at these temperatures, so-called warm-hot intergalactic 
medium (WHIM), mainly absorbs and emits photons in the 
ultraviolet (UV) and X-ray wavelength bands. 

Searching for the missing baryons has been conducted ex- 
tensively in the U V band. Indeed, besides the cool phase 
Lya absorbers (e.g., Weymann et al. 1998; P enton et al.l2000l 
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2004), the absorption lines of C m, C iv, N v, O vi, Si m, and 
Si iv have been routinely observed in spectra of background 
quasi-stellar objects ( QSOs; e.g., Itripp et al.l 12000. 2008; 



Prochaska et al 
2005L 120081:" 
Yao et all 120 HIT 



2004t lLehner et al.l|2006t iDanforth & Shull I 
Danforth et all 120061: iThom & Chen I 120081; 

that a significant amount of 
baryons exist in the highly ionized absorbers. These high 
ions, O vi in particular, are believed to trace the WHIM at 
the low-temperature end (~ 10 5 - 10 6 K). However, all these 
ions can also be produce d through photoionization in the in- 
terg alactic environme nt ( Oppenhei mer & Dav e 2009; but see 
also lTepper-Garcfa et al.l201ll and lSmith et al.l2012l for a con- 
trary view) and thus some of them may not contribute to 
the canonical WHIM (i.e., the shock-heated IGM expected 
from simulations). How many of these absorbers originate 
in the real WHIM is still under debate (IDanforth & Shull I 
120081; iTripp et al.ll2008l) . Ne vm and broad H i Lya absorbers 
(BLAs) are complementary and valuable tracers of hot ab- 
sorbing WHIM in UV bands. However, there a re only a few 
detec t ions with marginal significanc e of Ne vm (ISavage et al.1 
120051: iNarayanan etaD.120091 120111) and detections of BLAs 
require very high signal-to-noise-ratio (S/N) of s pectra to re- 
solve broad and weak signals from the continua dRichter et al] 
1 20061: iLehner et al.ll2007[ ISavage et al.1 1201 It IDanforth et all 
120101,1201 ll) . Converting measurement of these UV absorbers 
to census of baryonic matter also depends on ionization frac- 
tions and metallicity of the IGM, which in many absorbers are 
still poorly known. Nevertheless, BLAs and O vi-absorbers, 
although partially overlapped with photoionized Lya ab- 
sorbers, are estimated to contribute an additio nal ~ 25% to the 
baryon inventory at the current epoch (e.g.. IDanforth & Shull I 
200| IDanforth et aDl20TTt IShull et alJl201 ll) . 

X-ray observations of the intergalactic O vn and O vm 
emission/absorption features could provide essential informa- 
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tion for establishing the existence of the WHIM and com- 
pleting the baryon inventory in the local universe. At tem- 
peratures r>1 6 K, hydrogen and hel ium are nearly com- 
pletely ionized (Gnat & Sternberg 2007). Without metals, the 
thermal gas emits/absorbs photons through bremsstrahlung, 
which unfortunately i s optically thin even at intergalactic 
scales (Bregma n 120071) . Because the K-shell transitions of all 
elements heavier than lithium and the L-shell transitions of el- 
emen ts heavier than iron lie in X-ray bands (Paer els & Kahnl 
2003), metals can greatly increase the emissivity of the gas. 
However, because of the density-squared dependence, the 
WHIM emission is expected to be weak and difficult to detect. 
In fact, the IGM emission has been directly observed only 
from t he dense regions l i ke intracluster and intragroup media 
(e.g.. iWang et al.lll997t iMulchaev I l2000t lAllen et alJ l200l 
ISun et all |2009). The WHIM emission signals were also 
claimed in the angular correlation of the diffuse X-r ay back- 
ground (ISoItan et all 19991: lUrsino & Galeazzi 120061) . but dis- 
entangling the real WHIM signal from the unresolved point 
sources and loca l Galactic hot diffus e gas has been a major 
uncertainty (e.g., Galeazzi et al. 20071). 

Unlike emission lines, absorption lines measure the column 
densities and thus directly sample the total mass of the inter- 
vening gas along a line of sight. Oxygen is the most abundant 
metal element. In contrast of O vi as a minority ionization 
state, O vii is the most abundant ion for a gas at temperatures 
of T ~ 10 5 5 - 10 6 3 K and O vm takes over when T > 10 6 3 
K (Gnat & Sternberg 2007). Thus they have advantages over 
O vi in estimating the total baryon contained in highly ionized 
absorbers. O vn and O vm, whose production ionization po- 
tentials are 138.1 eV and 739.3 eV compared to 113.9 eV of 
O vi, are also hard to p roduce through photoionization in th e 
IGM environment (e.g. JCen & Fang l2006t IChen et alj2003l) . 
Theoretical calculations indicate that the O vn column density 
is about 10 times higher than O vi in a shock-heated gas (i.e., 
ypvii > KWqvi) wherea s Afovii < 3Afovi in a photoionized gas 
dFurlanetto et al. 1 120051) . Therefore observations of the IGM 
O vn and O vm absorption lines are crucial, not only for con- 
straining the properties of the WHIM at high temperature end 
but also for probing the nature of the commonly observed 
O vi-absorbers. 

Unfortunately, most attempts at searching the X-ray WHIM 
absorption lines have been frustrated. O vii and O vm ab- 
sorption lines at z = have been unambiguously detected in 
many background QSOs. However, all these X- ray absorp- 
tions are inconsistent with an intergalactic origin ( Fang et all 
l2006t iBregman & Llovd-Davies Il2007t lYao et alJEOlOb . but 
rather can be attributed to the Galactic diffuse hot gas (e.g., 
Yao & Wang 112001 12001: lYao et al.ll2008l 120091: [Wang et al. 



20051) . Perhaps the most compelling intergalactic result is 
the detection of O vii absorption at the redshift of the large 
galaxy structure Sculptor Wall (z a b s = 0.03) along the QSO 
sight line H 2356-309 (zgso = 0.165), albeit of low signifi- 
cance (~ 3a; iBuote et alj|2009t iFang et alj|20!ob . However, 
because of the high number density of galaxies along the sight 
line, the absorption may sample the halo gas of one or more 
intervening galaxies with small impact distances (<50 kpc; 
IWilliams et al.l l2010). i.e., mimicking the O vii and O vm ab- 
sorption lines at z = 0. Therefore, it may not be representa- 
tive of the typical WHIM. Furthermore, attempts at searching 
for the similar absorption features at the redshift of another 
large structure, Pisces-C etus along the same sight line, failed 
(Zappacost a et alj|2010h . All other claimed WHIM O vii and 



O vin absorptions have been highly debated. For instance, the 
detected O vm absorption line at z = 0.0554 in Ch andra spec- 
trum of PKS 2155-304 (IFang et al.ll2002aL 120071) cannot be 
confirmed by the XMM-Newton observations (Cagnoni et al. 
2004), although at nearby redshift a small galaxy group 
and the correspondi ng UV O vi and Lyq a bsorption lines 
have b een identified dShull etalj [l998. 2003). iNicastro et all 
(120051) obtained a spectrum with unprecedented S/N during 
the burst states of Mrk 421 and reported two WHIM detec- 
tions at z = 0.011 and z = 0.027. Again, observations with 
XMM-Newton cannot confirm the dete ctions an d reported 
significances have also been questione d ( Kaastra 



id reported 
etaDl2006; 



iRasmussen et alj2007l) . IWilliams etail d2010l) recently found 
a galaxy filament at z = 0.027, which makes the reported O vii 
WHIM detections along the Mrk 421 sight line back to be a 
subject of discussions. 

In this work, we aim to provide detection limits of current 
X-ray observatories, Chandra and XMM-Newton, in measur- 
ing the WHIM through X-ray absorption lines, and establish 
requirements (e.g., effective area and spectral resolution) for 
next-generation X-ray telescopes. We begin our investigation 
by scrutinizing the controversial WHIM detections along the 
Mrk 42 1 sight line. We extensively explore all the available 
Chandra and XMM-Newton observations, and then run boot- 
strap simulations to test the reliability of any absorption fea- 
ture. 

The paper is organized as follows. In Section [2] we de- 
scribes the observations and our data reduction processes, and 
in Section [3] we report the data analysis results and compare 
them with previous work. In SectionjU we run bootstrap sim- 
ulations to explore the detection limits of Chandra and XMM- 
Newton and to establish requirements of next-generation X- 
ray telescopes. We summarize our results in Section|5] 

Throughout the paper, we adopt the ato mic data from 
iVerner et al.l (fl99 6) and solar abundances from Asplu ndet al.l 
(2009). We conduct ou r data analysis within XSPEC (ver- 
sion 12.6: lArnaud I[l9 96), and report the la confidence range 
or 3a limit for a fitting parameter. We take the flux nor- 
malization of 1 Crab = 2 x 10" 8 erg s" 1 cm" 2 at 0.5-2.0 keV. 
Thus 0.2 mCrab is 2.15 x 10" 13 erg s _1 cm -2 A -1 for a flat 
spectrum in wavelength space, which corresponds to 2.34 x 
10" 4 photons s" 1 cm" 2 A" 1 around the O vn Ka (21.602 A) 
at the rest-frame. We define measurement significance level 
(SL) of an absorption line as 



SL = 



EW 
AEW' 



(1) 



where EW and AEW are the equivalent width and its la un- 
certainty for an absorption line. 

2. OBSERVATIONS AND DATA REDUCTION 

We used all grating observations of Mrk 421 observed with 
both Chandra and XMM-Newton that are available as of 201 1 
April. Since Mrk 421 is a calibration source for both obser- 
vatories, observations were taken on a nearly regular basis. 
The Chandra data were downloaded from the CXO archive 
web siteQ, which included 26 observations taken with the low 
energy transmission grating (LETG; ObsIDs of 1715, 4148, 
4149, 5171, 5318, 5331, 5332, 6925, 8378, 8396, 10664, 
10665, 10667, 10668, 10669, 10671, 11605, 11606, 11607, 
11966, 11970, 11972, 11974, 12121, 12122, and 13097) with 

3 http://cda.harvard.edu/chaser/mainEntry.do 
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total exposure of 752.8 ks. This dataset does not include ob- 
servations taken with imaging mode and those with the high 
energy transmission grating (HETG). A quick estimate indi- 
cates that these excluded observations would contribute addi- 
tional < 10% to the total spectral counts, and thus our results 
will not be noticeably affected with or without them. For each 
employed observation, we used software package CIAO (ver- 
sion 4. 3)0 and followed the standard procedure to re-calibrate, 
extract spectrum, and calculate the instrumental redistribution 
matrix file (RMF) and ancillary response file (ARF). For ob- 
servations that used the Advanced CCD Imaging Spectrom- 
eter (ACIS) as the detector, we only utilized the first grating 
orders. The High Resolution Camera (HRC) does not have 
spectral resolution itself, and the spectral orders in observa- 
tions with HRC cannot be sorted out. Therefore, for these 
observations, we calculated the RMFs and ARFs up to six 
grating orders and combined them as one in spectral analy- 
sis to account for contributions from different orders (see de- 
scriptions in Section^. To imp rove spectral S/N, w e utilized 
the user-developed IDL tools (lYao & W ang 2007|) to com- 
bine the negative and positive spectral orders, coadd spectra 
of different observations, and calculate the averaged RMF and 
ARF by using the spectral counts around 20 A as weights. 
The claimed two WHIM detec tions were based on a subset of 
the observations listed above (Nic astro et al J 120051 : e.g., Ob- 
sIDs of 1715, 4148, and 414®. To scrutinize the reported 
WHIM detections, we produced three coadded Chandra spec- 
tra and the corresponding RMFs and ARFs, based upon ob- 
servations used by Nicastro et al. (Spectrum I), observations 
not included in their work (Spectrum II), and all observations 
together (Spectrum III), respectively. Figures QJ3] show six 
portions of these spectra. 

The total 76 XMM-Newton observations were downloaded 
from the XMM-Newton science archive Q. We used script 
package SAS (version 10.0)Qto reprocess these observations. 
To avoid the severe contamination from the background flares, 
for each observation we calculated the count rate on the chip 
CCD9, and created the good time intervals (GTIs) by exclud- 
ing the time periods with rates higher than 0.2 counts per sec- 
ond. After this filter, there are 54 remaining observations, 
and each of them has the integrated GTIs longer than 5 ks. 
We then extracted spectra, calculated RMFs for each of ex- 
posures taken with the reflection grating spectrometer (RGS), 
and combined the spectra of these 54 observations and the 
corresponding RMFs by running the script rgscombine to get 
the final coadded spectrum (Spectrum IV; Figure 0), which 
has a total exposure of 1.2 Ms. 

We used a binsize of 25 mA to re-bin all the coadded spec- 
tra, which is equivalent to an oversample factor of 2 for a 
^50 mA resolution element of the LETG and RGS. Table Q] 
summarizes the S/Ns of these spectra. 

3. SPECTRAL ANALYSES, RESULTS, AND DISCUSSION 

Our goal in this section is to exam ine the WHIM detec tions 
at redshifts z = 0.01 1 and z = 0.027. iNicastro et all d2005l) re- 
ported the WHIM absorption from transitions of Ne ix Ka, 
O vn K/3, O vm Ka, O vn K/3, N yii Ka and C iv Ka (Ta- 
ble 2 and Figure 8 in [Nicastro et al.ll2005l), whose rest-frame 

4 http://cxc.harvard.edu/ciao/index.html 

5 Nicastro et al. 12005) in fact also included two additional HETG obser- 
vations, which contributed negligible spectral counts to their final spectrum. 

6 http://xmm.esac.esa.int/xsa/index.shtml. 

7 http://xmm.esac.esa.int/sas/. 



wavelengths are 13.44 7, 18.629, 18.969,2 1.602, 28.787, and 
33.736 A, respectively (iVerner et alj|1996l) . Therefore, we fo- 
cus our attention on spectral ranges covering these transitions. 

Obtaining a good continuum model is crucial to absorption- 
line measurement. Since all coadded Chandra spectra con- 
tain the contribution from HRC observations, to account for 
grating order-overlapping issues (Section [2]), we decide to fit 
a broad range of the spectra from 7 to 40 A, which cov- 
ers the transitions to be examined. For a ease of compari- 
son of three Chandra spectra, we jointly fit the Spectra I, II, 
and III. We first use a Galactic-absorption modified power- 
law (wabs*powerlaw in xspec) to fit the continuum emis- 
sion, with neutral hydrogen column density Nh linked to- 
gether while the power-law index T and normalization are 
allowed to vary. We obtain an unacceptable fit with y 2 
1 1 646 over 3590 degrees of freedom (DOF), and we find 
that some "broad" features cannot be accounted for by this 
simple model. These features are mainly due to imperfect 
instrumental calibration around node boundaries, CCD chip 
gaps plus dithering effects, the oxygen absorption edge, and 
imperfe ct cross calibration between ACIS and H RC obser- 
vations dMarshall et al.ll2004l: INicastro et alJUooBI) . We then 
use Gaussian profiles to compensate the uncounted calibra- 
tion residuals. To minimize the effect of the known strong 
ISM absorption lines on continuum modeling, we include nar- 
row absorption lines (with widths a < 50 km s" 1 ) of O i, O n, 
C vi, O vi, O vn, and Ne ix Ka and O vn K/j transitions at 
their rest-frame wavel engths (IVerner et al.l Tl996: Jue tt et all 
l2^06tlYaoet"aIll2009h in our spectral fitting. The centroid 
energies and widths of these Gaussians are linked together 
while the normalizations are allowed to vary in the jointly 
spectral fitting. We find that, besides these narrow absorp- 
tion lines produced in the ISM, we need 10 broad and one 
narrow Gaussian profiles to obtain an acceptable fit, and lo- 
cations of these Gaussian are 10.53, 13.80, 14.60, 17.90, 
18.19, 19.02, 19.10, 23.24, 23.31, 23.66, and 29.89 A with 
widths (a) of 0.61, 2.56, 0.04, 0.00054, 0.17, 1.36, 0.10, 
1.45, 1.69, 0.41, and 0.39 A. The final X 2 = 4448.2 with 3512 
DOF and the constrained N R = 1 .092 ± 0.06 x 10 20 cm -2 , and 
T = 2.131 ±0.001, 2.011 ±0.001, 2.067 ±0.001 for Spectra 
I, II, and III, respectively. Because the absorption line mea- 
surement depends only on the local continuum, and because 
the centroids of these Gaussian profiles do not directly su- 
perpose on the wavelengths of the putative WHIM absorption 
lines, these Gaussian profiles will not have any effect on the 
line equivalent width measurement conducted below. XMM- 
Newton spectrum (Spectrum IV) does not have the grating- 
order-overlapping issue, so we use the wabs*powerlaw model 
to fit the six narrow spectral ranges as plotted in Figure [4] 
We exclude the bad pixels in our spectral analysis, and obtain 
x 2 /DOF = 1 .36 , 1 .49 , 1 .49 , 1 .28 , 1 .75 , 1 .56 with DOF of 1 37, 
75, 103, 105, 116 and 104, respectively. Figures [Q|4] show 
these best fit continua with the ISM absorption lines removed 
from spectral models. 

Now, let us examine the existence of the reported WHIM 
absorption lines. Figures [T]|4] show the same portions of 
Chandra and X MM-Newton spectra as shown in Figure 8 of 
INicastro et al.l (120051) . in which they demonstrated their detec- 
tions of the WHIM absorption lines at z = 0.01 1 and z = 0.027 
from transitions of Ne ix Ka, O vn K/3, O vm Ka, O vn Kf3, 
N vn Ka, and C iv Ka. In particular, Figure Q] reproduces 
their Figure 8, as they were extracted from the same Chan- 
dra observations. Several bad pixels of the RGS accidentally 
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FIG. 1. — Six portions of coadded C/iandra-LETG spectrum (green histogram) of Mrk 421 with the best-fit continuum model (blue curves). The model is 
co nvolved with the instr umental line spread functions. The binsize is 25 mA. The spectrum (Spectrum I) was extracted from the same observations as used 
by Nicastro et al. 12005), and the plots reproduce their Figure 8. Ka and K/3 transitions of eight ions covered by these spectral ranges are labeled in the plots, 
and the three vertical bars from left to right mark the corresponding rest-frame wavelengths (thin line), the expected wavelengths at z = 0.011 (medium line), 
and i — 0.027 (thick line). The red curve in panel (c) indicates the additional component added to local continuum in order to "amplify" the significance of the 
putative O VII Ka absorption line at z = 0.01 1. The additional component is statistically unnecessary, and the red curve is for comparion purpose only. See text 
in Section[3]for details. 



TABLE 1 
Spectral analysis results 













Chandra 






XMM-Newton 


line 


redshift 


S/N 1 


EW (mA) 1 


S/N 2 


EW (mA) 2 


S/N 3 


EW (mA) 3 


S/N 


EW (mA) 


O VII 4 


z = 


68.9 


10.9 ±0.9 


64.9 


12.5 ±0.9 


94.5 


11.7±0.6 


193.0 


14.9 ±0.6 




2 = 0.011 


67.6 


<3.0 


62.8 


< 1.5 


92.2 


< 1.6 








Z = 0.027 6 


66.2 


<3.1 


60.6 


< 1.5 


89.6 


< 2.0 


187.4 


< 1.6 


O VII 5 


Z = 


68.9 


11.3 ±1.0 


64.9 


11.8± 1.1 


94.5 


11.6±0.9 








z = 0.011 


67.6 


3.3 ± 1.2 


62.8 


< 2.0 


92.2 


2.0± 1.1 








z = 0.027 6 


66.2 


<3.3 


60.6 


< 1.5 


89.6 


< 2.0 







NOTE. — The S/N is given per 50 mA resolution element around the lines, and EW is the measured equivalent width. Error ranges are quoted at Icr significance 
level, while the upper limits are at 3<r. 

1 Measurements are obtained from the spectrum corresponding to observations used by Nicastro et al. 1 2005) (Spectrum I; Figure[T}. 

2 Measurements are obtained from the spectrum corresponding to those observations not used by Nicastro et al. (Spectrum II; Figure|2). 

3 Measurements are obtained from the spectrum corresponding to all observations (Spectrum III; Figure[3J. 

4 EW measurements in these three rows are obtained from the best-fit continua without compensation of the 1 1th broad Gaussian (blue curves in Figures [TBI . 

5 EW measurements in these three rows are obtained from the best-fit continua with compensation of the 11th broad Gaussian (red curves in Figures [7131 . 
However, since the 1 1th broad Gaussian profile is statistically unnecessary, these measurements are for comparison purpose only. 

6 The measured upper limits to EW of the O VII Ka are obtained by fixing the line centroid at z. = 0.028 instead of z = 0.027, since Nicastro et al. 1 2005) reported 
O VII line at the former redshift. 

See text for details. 
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FIG. 3. — Same as Figure[T] except that the spectrum (Spectrum III) was extracted from all available Oiandra-LETG observations. 



6 



Yao et al. 



3.0xl0 4 



2.8xl0 4 



2.6xl0 4 



1 


Ne XxS^' - 


- Ne X Kc, 


II 












1.4xl0 4 












N 




(e) : 


1.2x10* 






1 1 


l.OxlO 4 










8.0xl0 3 










6.0xl0 3 






I J 






28.5 29.0 29.5 33.5 34.0 34.5 

Wavelength (A) 

FIG. 4. — Same as Figure[T] except that the spectrum (Spectrum IV) was extracted from all available XMM-Newton-RGS observations. Channels/bins colored 
as grey were contaminated by hot pixels and were ignored in our spectral analysis. Panels (b)-(f) were extracted from the RGS1 observations, and panel (a) was 
extracted from the RGS2 observations. 



fall around the 21.8 A region and cause an instrument feature 
(Figure 0J. Although a careful calibration could yield a cor- 
rect instrumental response of the RGS, in this work we do not 
use the XMM-Newton spectrum to assess the WHIM detec- 
tion at z = 0.011. Visual inspec tion reveals that none of the 
above WHIM lines reported by iNicastro et all ((2005) is de- 
tected in the previously Chandra observations they analyzed 
(Figure Q~|l, and in the newly available Chandra observations 
(Figure|2), and in the XMM-Newton observations (Figure|4]l. 

The reported WHIM detections could be due to the im- 
proper continuum placement on the spectrum. Visual inspec- 
tion also reveals that there might be an absorption feature at 
~ 21.8 A in Spectrum I, which corresponds to the reported 
O vn Ka WHIM absorption line at z = 0.011 (panel (c) in 
Figure [TJ. However, such an absorption feature is not visible 
in Spectra II and III (panel (c) in Figures [2] and 0, and its 
significance in Spectrum I depends largely on how the local 
continuum is placed. Please note that our spectral profiles of 
the local continua were obtained by jointly fitting the Spec- 
tra I, II, and III (see above), and the best-fit model seems to 
account for the general spectral variation reasonably well. To 
quantitatively assess the existence of the O vn Ka WHIM line 
at z = 0.01 1, in the joint analysis of the Spectra I, II, and III, in 
addition to the 10 Gaussian profiles, we add the 1 1th "broad" 
Gaussian to the local continuum and another narrower Gaus- 
sian with width a = 50 km s" 1 to represent the WHIM ab- 
sorption line. Again, we link the line centroids and widths 
together in our joint fitting but allow the normalizations to 



Different widths yield nearly identical resul ts as long as the line is still 
unresolved by the instrument. Please see Section |4~Tl for justifications for the 
line width. 



vary among the Spectra I, II, and III. The best fit yields a line 
centroid of 21.90 A and a = 0.050 A for the "broad" compo- 
nent and a centroid of 2 1 . 87 A for the narrow component. The 
red curves in panel (c) of Figures[T][3]show the modified local 
continua. These additional components reduce the \ 2 of 26.4 
by adding nine DOF in total. However, applying the best-fit 
models to individual spectra yields a \ 2 change of 13.7, 4.3, 
and 10.4 by introducing five additional DOF to spectral fit- 
ting of Spectrum I, II, and III, respectively. All these Chandra 
spectra can be well described by the same continuum profile 
with various normalizations, which has been proved to be rea- 
sonable to other wavelength ranges (panels (a)-(b) and (d)-(f) 
in Figures [HO. Because we find essentially no improvement 
in fitting the Spectrum II, adding the "broad" Gaussian profile 
to the local continuum around 21.9 A is statistically unnes- 
sary. 

We next attempt to measure the equivalent widths (EW) of 
the putative WHIM absorption lines. Since we have not con- 
sistently detected any WHIM a bsorption of any transitions re- 
ported by INicastro et alJ (120051) . we focus here on measuring 
the upper limit to the equivalent width (EW) of the O vn Ka 
line, which is expected to be the most abundant transition in 
a broad temperature range of the WHIM (SectionQ]). At red- 
shifts z = 0.01 1 and z = 0.027, the rest-frame O vn Ka at 21.60 
A is shifted to 21.84 and 22.19 A, respectively, and therefore 
our measurements are obtained from the spectral range de- 
scribed by the continuum colored as blue in panel (c) of Fig- 
ures Q]|4] T o make a fair comparison with the results obtained 
by Nicastr o~et al.l d2005l) . we also make similar measurements 
by using the continuum colored as red. We should empha- 
size that the latter measurements are only for comparison pur- 
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poses, since the additional component to the local continuum 
that "amplifies" the significance of the absorption line at 21.8 
A is statistically unnecessary. Our results are reported in Ta- 
ble!]] 

In summary, our analysis and results do not support the ex- 
istence of the two WHIM systems at z = 0.01 1 and z = 0.027. 
We analyze all the available Chandra and XMM-Newton ob- 
servations of Mrk 421, and have not seen any consistent 
WHIM absorption of any transition at eith er z = 0.011 or 
z = 0.027 reported by Nicas troet al.l ([2005 ). We measure 
the EW of the O vn Ka absorption line. For the system at 
z = 0.027, we obtain a firm 3a upper limit of EW < 1.5 mA 
from the spectrum (Spectrum II) extracted from the newly 
available Chandra observations, in contrast to EW = 2.2 ±0.8 
mA obtained by Nicastro et al. For the system at z = 0.01 1, 
we confirm the existence of a s mall dip in the spectrum (Spec- 
trum I) identical to that used by Nicastr o et al.l ( 120051) . but the 
similar spectral feature is absent in Spectrum II. The measure- 
ment of the line heavily depends on how the local continuum 
is placed, and the additional spectral component that amplifies 
its significance is statistically unnecessary. We obtain an up- 
per limit of EW < 1 .5 mA from the new spectrum, in contrast 
to EW = 3.0^1 mA reported by INicastro et al.l (120051) . The 
findings in this work are consi s tent w ith the previous inves- 
tigation made by Kaastra et al.l d2006l) . in which the authors 
scrutinized the same data set used by Nicastro et al. and a sub- 
set of calibration observations of Chandra and XMM-Newton. 
Therefore, we conclude that there is no WHIM line detected 
at either z = 0.01 1 or z = 0.027 along the Mrk 421 sight line. 

There are also other important absorption lines in the spec- 
tral range of interest (Figures HE). The prominent lines are 
O vn, Ne ix, and C iv Ka and O vn K/3 at z ~ 0, and the O vi 
Ka with rest-frame wavelength 22.0 A at z ~ is also visi- 
ble in all four spectra albeit less significant. Again, we only 
measure the EWs of O vn Ka and report them in TableQ] The 
O vn line at z ~ can be well explained as abso rption from 
the G alactic diffuse interstellar medium (ISM) (lYao & Wang I 
120071), and the O vi line i s also believed to originate from the 
ISM (ISavage et al.ll2005l) . It is worth noting that the EW of 
O vn Ka at z ~ measured from XMM-Newton observations 
is > 3c larger^ than that obtained from Chandra observations 
(Table [0. Such a discrepan cy in O vn has also been mea- 
sured by Kaast ra et al.1 (1200 6). The causes to this discrepancy 
are still under investigation, and the results will be published 
elsewhere. In this paper, we mainly focus on X-ray absorp- 
tion measurements of the WHIM and will not discuss these 
non-WHIM measurements any further. 

4. BOOTSTRAP SIMULATIONS 

In Sections [2] and [3] we reprocessed all the available obser- 
vations and obtained two spectra with the highest S/Ns ever 
collected with Chandra and XMM-Newton. Despite the in- 
creased S/N, we s till cannot confirm the X-ray WHIM detec- 
tions reported by INicastro et al.l (120051) . Given such a frus- 
trating situation (see also discussions in Section [T), several 
crucial questions need to be answered before further major 
time allocations search for the WHIM through X-ray absorp- 
tion line measurement: How high S/Ns are needed to firmly 

9 The differential significance is calculated as 

\EWi-EW 2 \/V(AEW [ ) 2 + (AEW 2 ) 2 , 
where AEW is the 1<t error of EW measurement. 



establish a WHIM detection? What is the O vn column den- 
sity range suitable for Chandra and XMM-Newton to explore 
with reasonably long exposure time? Is it possible to conduct 
a systematic study of the WHIM using the current X-ray tele- 
scopes? If not, what are the requirements for next-generation 
X-ray telescopes? 

4. 1 . Expected WHIM properties and required instrumental 
sensitivities 

Before running bootstrap simulations to address these ques- 
tions, we first review the WHIM properties based on the 
results of numerical simulations. We then estimate the re- 
quired instrumental sensitivity to detect the WHIM by an- 
swering the following two specific questions: First, what 
O vn column density is needed to probe the majority of the 
baryons contained in the O vn-bearing WHIM? Cosmolog- 
ical simulations predict different O vn column density dis- 
tributions per unit redshift with different r ecipes for galactic 
feedback and intergalactic photoioniza tion dFang et al 12 002b: 
ICen & Fang Il2006[ ISmith et al.ll2012h . In Figure [5] we con- 
vert two recent distributions (Figure 5 in lCen&Fang| |2006 
and figures in Smith et al] 120121) to the fraction of baryons 
contained in the O vn-bearing WHIM as a function of O vn 
column density. In this conversion, we assume a single 
temperature and uni formly distr i buted metallicity for all ab- 
sorbers. However, Smi th et al] ((2012) find covariance be- 
tween Afovii and metallicity, and the WHIM is expected to 
be multiphase. To obtain a more realistic distribution, one 
should consider the ionization and metallicity correction to- 
gether. However, without a detailed description of the WHIM 
properties (which in fact are our final goals), such a correction 
is difficult to apply. Considering the general trend that higher 
column density absorbers have higher metallicity, the plot- 
ted baryon fractions contained in the higher column absorbers 
(e.g., A^ovn^lO 15 cm" 2 ) should be regarded as upper limits. 
Furthermore, since th ere is no clear sign of convergence a t 
low column densities (ICen & Fang Il2006t ISmith et al.ll2012h . 
in our conversion we use a 4- or 6-degree polynomial to 
fit the number density distributions, which artificially gener- 
ate turnover at column densities 10 10 - 10 12 cm" 2 . The low 
column density absorbers could contain even more baryons 
if the distributions converge at lower columns. Neverthe- 
less, Figure [5] gives an estimate of the baryon fraction dis- 
tribution, which indicates that ^55% (or 65% or 80%, de- 
pending on which simulation is selected) of baryons are con- 
tained in absorbers with Afovn > 10 14 5 cm" 2 and that ^30% 
(or 45% or 55%) of baryons are contained in absorbers with 
JVovn > 10 15 cm" 2 . 

Another important, yet to be determined property is the dis- 
persion/Doppler velocity b v of the O vn-bearing gas. This 
property, although not essential f or th e detectability of Chan- 
dra and XMM-Newton (Section l4~2l . is important for esti- 
mating the req uired spectral resolution of future X-ray instru- 
ments (Section PO) . Here we use O vi as a reference. Assum- 
ing most detected IGM O vi tracing the WHIM at tempera- 
tures ~ 10 5 5 K and the total b v being a combination of ther- 
mal and turbulent motions through quadrature b\ = b^ h + bf urb , 
the median b v ~ 30 km s" 1 of the O v i absorption lines 
dDanforth & Shull 1120081: iTripp et al.ll2008l) suggests a turbu- 
lent velocity of Z? t urb ~ 24 km s" 1 for the O vi-bearing gas. 
Higher temperature gas may have higher turbulent velocity. 
Thus b v of the O vn-bearing gas at temperatures ~ 10 6 K is 
expected to be > 40 km s" 1 . Another useful reference is the 
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FIG . 5 . — Fraction of baryons contained in the O VH-bearing gas as a func- 
tion of O VI I colum n density. The solid line is derived from Figure 5 in 
Cen & Fang (2006). The dashed and dotted lines are derived from figures in 
Smith et al. (2012), which consider collisional ionization (CIE) only and CIE 
plus photoionization, respectively. See text for details. 

hot phase Galactic ISM traced by the O vn absorption line at 
z ~ 0, although the IGM environment could be substantially 
different from the ISM. Nevertheless, along the Mrk 421 and 
LMC X-3 sight lines, the b v of the O vn gas can be as large as 
100 k m s" 1 dWang et al.H200l lYao & Wang 1120071 lYao et all 
2009). We therefore will explore three b v values, 50, 75, and 
100 km s" 1 , in our simulations. 

Second, what is the flux limit that X-ray telescopes should 
reach in order to statistically investigate the WHIM? The an- 
swer to this question depends largely on how many WHIM 
systems are planned to detect. Clearly when more systems 
are detected, one can investigate WHIM properties, like ther- 
mal and dynamic properties, metallicities, and their evolution 
with redshift. To obtain a spectrum with a requested S/N, the 
required exposure time is inversely proportional to the flux of 
a QSO. And the higher the redshift of a QSO, the higher pos- 
sibility there is intervening WHIM absorbers along the line 
of sight. Therefore, to estimate the number of WHIM sys- 
tems that could be detected, one should create a target list by 
considering the QSO redshift and flux as a whole. 

In this work, instead of providing a target list for the 
WHIM investigation, we require that future X-ray tele- 
scopes be able to observe and search for the corresponding 
O vn absorbers along the majority of the sight lines where 
IGM O vi absorbers have already been reported (Sec- 
tion [p- In this regard , we examine the targ et lists in 
lDanforth& Shulll (120081) and iTripp et alj (120081) and find 
that there are ~ 80 IGM O vi systems along 28 QSO sight 
lines reported so far. Among them, we have assembled the 
available soft X-ray fluxes for 24 sight lines from the 
bright source catalogs of Einstein, ASCA, BeppoSAX, and 
ROSAT, and from severa l pointing observations with X MM- 
Newt on and Chandra dWilkes & Elvis I 119871 lUedaetaTl 
20051: lReeves& Turner! I2000t IVerrecchia et al.1 120071: 
Bianc hi et al.l 120091 iPaoadakis et al.1 120071; iFoschini et alJ 
2006] iPorquet et alJ 120041 IWhite et all l2000t iLeighly et all 
2007). These fluxes range from 0.006 to 46.5 mCrab with a 
median of ~ 0.3 mCrab. Fifteen of them have fluxes >0.2 
mCrab and sample a total redshift pathlength of Az = 2.1. 
For the WHIM O vn -absorbers with JVovu > 10 15 cm" 2 , 
ICen & Ostriker I d2006|) predicted t he ab sorber frequency to 
be dM/dz ~ 7, and ISmith et all (120121) predicted the fre- 

10 We use a power-law form f(E) ~ E~ l A dE to convert the reported flux 
to the energy range of 0.5-2.0 keV (Section[T}. 
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FIG. 6. — Curve of growth of the equivalent width of O VII Kq absorp- 
tion line versus O VII column density with three different Doppler dispersion 
velocities (b v ). 

quency to be ~ 20 considering collisional ionization only and 
dM I dz ~ 5 if photoionization is included. For absorbers with 
Afovn > 10 14 5 cm" 2 , dM/dz is predicted to be ~ 18, 12, and 
110, respectively. Therefore, X-ray telescopes should be able 
to facilitate absorption line study from background sources 
with flux of ~ 0.2 mCrab, and >10 WHIM systems with 
Novn > 10 15 cm" 2 and >24 systems with JVovu > 10 14 ' 5 cm" 2 
are expected in these sources along which the corresponding 
O vi absorptions are already known. 

4.2. Detection limits of current X-ray telescopes 

We now use bootstrap simulations to examine the de- 
tectability of Chandra-LETG and XMM-Newton-RGS. In our 
simulations, w e utilize the user-d eveloped absorption line 
model absline (lYao & Wang 11200*51) to model O vn absorbers 
with column density jVovii an d dispersion velocity b v , and use 
a PL continuum to model an emission spectrum. With the 
input parameters of Novu, b v , PL normalization, and an ex- 
posure time, we simulate two spectra within Xspec, using the 
RMFs and ARFs produced in Section [21 to mimic Chandra 
observations and XMM-Newton observations, respectively. 
From the simulated spectra, we use a Gaussian profile to rep- 
resent the absorption line and measure its EW and uncertainty 
AEW at the expected wavelength position and then calculate 
the line SL (Eq. [2]). For a set of parameters, we repeat this 
procedure 1000 times and record the measurement SL of the 
line. We then sort the SL and select the 100th value (i.e., 
90% simulation trials yield the SL equal to or higher than the 
selected value). We make the same measurement for various 
combinations of different exposure times, b v , and A^ovn- Since 
the full-width-half-maximum (FWHM) of the LSF of both 
Chandra-LETG and XMM-Newton-RGS is -50 mA, which 
is much broader than the expected b v values of the WHIM 
O vn absorbers (Section 14.11 ). all simulated absorption lines 
appear unresolved. What really matters for measurement is 
the effective EW of an absorption line. To provide a direct 
comparison to observations, we convert the combination of 
Afovn and b v to EW through the curve-of-growth (COG), and 
the combination of exposure time, RMF, and ARF to S/N per 
50-mA resolution element. For reference, we plot the COG 
of O vn Ka line in Figure [6] 

Figure [7] shows our simulation results. Taking the simula- 
tions for Oiandra-LETG for an example, our results can be 
interpreted as: from a spectrum with S/N = 50 per resolution 
element, an O vn Ka absorption line with EW = 3.4, 5.7, 8.1, 
10.3 mA can be measured at >2, 4, 6 and 8er SLs, respec- 
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tively, in 900 out of 1000 observations. In another word, to 
measure the absorption line with these EWs at 2, 4, 6 and 8cr 
or higher SLs, a spectrum with S/N > 50 is required. Interest- 
ingly, the measured line significances behave nicely like PL 
functions in the EW-S/N space. For reference, we formulize 
them as, 

3 

log(W) = a x log(S/A0 + Pi x SV, (2) 

i=0 

where SL is 2, 3, 4, 5, 6, 7, or 8. For Chandra simula- 
tions, a = -0.932, = 1-79108, 0.188439, -0.0148031, 
and 0.000476159, and for XMM-Newton, a = -0.970, /3 -3 = 
2.14183, 0.217703, -0.0197475, and 0.000748764. 

Let us use our simulation results to further examine the two 
putative W HIM systems. I n Figure [7] we also plot the EWs 
reported by iNicastro et al.l (12005b and their 1 a errors at the 
corresponding S/N positions of Spectrum I. We find that their 
reported SLs appear to be consistent with expectations of our 
simulations. However, we stress again that their SL measure- 
ment could have been largely impacted by the improper con- 
tinuum placement (Section Our simulations indicate that 
an O vii Kg absorpt ion line with EW=3.0 mA at z= 0.011 
(INicastro et al.l l2005) should be measured at > 3.7er SL from 
Chandra Spectrum III with S/N- 92 (Eq.E Figure^). This is 
in contrast to the fact that we only obtain an upper limit (Fig- 
ured] TableQ]). At z = 0.027, an O vn line with EW = 2.2 mA 
should be detected at >2.3er from the coadded Chandra and 
XMM-Newton spectra, which, again, is in contrast to the fact 
that we only obtain upper limits. We therefore conclude that 
the putative O vn WHIM absorptions, if they exist at all, must 
have lower EWs than reported. 

Is the reported WHI M absorption at z = 0.0 27 associated 
with a galaxy filament? Williams et al.l (120101) discovered a 
filament of late-type galaxies at z = 0.027, suggesting a pos- 
sible association between the reported WHIM absorber and 
the filament. In the spectrum of PKS 0405-123 observed with 
the Cosmic Origins Spectrograph (COS) aboard the Hubble 
Space Telescope, a similar association has also been sug - 
gested by a broad O vi absorption line (Savage et al.H2010l) . 
The authors estimated that the WHIM contained within the 
abundant late-type galaxy groups could make up to 15% of 
the total baryons in the local universe. However, the non- 
detection of O vii at z = 0.027 in this work leads us to con- 
clude that there is, if at all, much less O vii than reported. A 
recent search in a high-S/N COS spectrum of Mrk 421 also 
failed to find any ab sorption feature of BLA at the redshift 
(iDanforth et alJl20llT) . Similarly, there is no dete cted hot in- 
tragroup O vii-bearing gas in our Local G roup dFang et alj 
l2006t lYao& Wang I [20071 iBregmanl 120071) . Therefore, if 
there is any gas associated with the galaxy filament at z = 
0.027 along the Mrk 421 sight line, it must be at higher tem- 
peratures (T > 10 6 5 K) and/or be metal poor and thus contain 
too little O vii and H i to be detected. 

Our simulation results also have important implications for 
using current X-ray observatories to search for the WHIM 
absorption. First, Chandra and XMM-Newton can probe 
the WHIM absorbers with Afovn^lO 16 cm" 2 with exposure 
times of — 1 Ms from a QSO spectrum with a flux of — 0.2 
mCrab. The O vii column densities correspond to EWs of 
> 10 — 15 mA (Figure [6j, which require a Chandra (XMM- 
Newton) spectrum with S/N = 18-27 (41-63) to measure 
it at >4cr significance levels (Figure [7J Eq. [2). To col- 



lect a spectrum with such S/N, takes 1.4-3.1 Ms for Chan- 
dra (adopting A = 20 cm 2 effective area) and 3.6-8.5 Ms for 
XMM-Newton (adopting A = 40 cm 2 effective area) from a 
QSO with 0.2 mCrab flux. This is consistent with the com- 
monly detected O vii absorption lines, which trace the Galac- 
tic hot ISM, in the spectra of l ocal QSOs (z < . 1) with fluxes 
> 2 m Crab dFang et al.ll2006t iBregman 112001 lYao & Wang] 
120071) . However, the simulated dAf/dz of O vii drops by 
about two orders of magnitude between iVovn > 10 15 cm -2 
and Afovii > 10 16 cm" 2 dSmith et al.ll2012l) . meaning that such 
high column O vn-absorbers are extremely rare and may only 
exist in the circumgalactic e nvironment around galaxy struc- 
tures (e.g. JFang et alJl20ToT) . 

Second, it is likely infeasible to conduct a systematic study 
of the WHIM with Chandra or XMM-Newton. The EW for 
Mdvii = 10 15 cm -2 is — 2.5 mA (Figure [6J, which requires a 
spectrum with S/N- 120 (Chandra) and 260 (XMM-Newton) 
to measure it at — 4<r significance level (Figure [7). It would 
require 61 Ms for Chandra and 144 Ms for XMM-Newton to 
collect a spectrum with such a high S/N from a QSO with 0.2 
mCrab flux. Even for a blazar at its burst state like Mrk 421, 
it still takes 0.6 (1.4) Ms Chandra (XMM-Newton) time from 
the source continuously being bright at 50 mCrab level. It 
would take 341 Ms for Chandra and 804 Ms for XMM- 
Newton to complete the survey of the 15 background QSOs 
with fluxes >0. 2 m Crab along which IGM O vi has been re- 
ported (Section l4~Tl i. Furthermore, since the expected WHIM 
O vii absorption lines cannot be resolved, it is impossible to 
study the dynamics of the WHIM with either Giandra-LETG 
or XMM-Newton-RGS. Clearly, we need new X-ray spec- 
trographs with much improved spectral resolution and much 
larger collecting area to systematically study the WHIM. 

4.3. Requirement for next generation X-ray telescopes 

The goal in this section is to establish the required spec- 
tral resolution and effective area for next generation X-ray 
telescopes. Since there is no existing X-ray telescope with 
spectral resolution better than those of Chandra and XMM- 
Newton, we use Gaussian profiles to simulate the LSFs of the 
putative X-ray spectrographs. We assume a spectral resolu- 
tion, which is defined Q as R = A/ A A, of 1500, 3000, 4000, 
and 6000, respectively, with effective area of A = 1000 cm 2 
throughout the whole wavelength coverage. We use four spec- 
tral bins within a FWHM spectral range (i.e., an oversample 
factor of 4) in our simulations. With a set of selected R, Afovii, 
and b v of the O vn-beari ng ga s, we follow the same procedure 
as described in Section 14.21 to simulate an absorption spec- 
trum of an background source with 0.2 mCrab flux for an as- 
sumed exposure time, and then measure the EW and AEW 
from the simulated spectrum. We repeat the procedure 1000 
times and select the SL value so that 90% trials have the same 
or higher SL. In a real observation, background source flux, 
instrumental effective area, and exposure time can compen- 
sate each other, so we define their product (FAE) as a figure 
of merit. Figure^shows our simulation results. For reference, 
we also label the corresponding EW to the labeled Afovn at top 
X-axis of the plots. Table[2]summarizes the required exposure 
time for the future spectrograph with 1000 cm 2 effective area 
and different spectral resolutions to constrain the O vn ab- 
sorption with Afovii = 10 15 cm -2 and three different b v values 
at > 4 and 6cr significance levels from a background source 

1 1 We use the FWHM of the LSF as AA. 
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Dispersion velocity b v 




R 


50 km s" 1 


75kms-' 


100 km s" 1 


1500 


416.9 (808.8) 


398.8 (794.6) 


414.6 (799.5) 


3000 


217.4 (433.9) 


237.6 (492.6) 


277.8 (549.9) 


4000 


171.1 (350.3) 


212.1 (425.7) 


241.2 (523.0) 


6000 


126.3 (276.5) 


162.9 (368.2) 


208.3 (461.7) 



FIG. 7. — The 2, 4, 6, and 8a detection significance of O VII Ka versus S/N per 50-mA resolution element with Oiandra-LETG (a) and XMM-Newton-RGS 
(b). The vertical red and blue lines mark the S/N levels of the final coadded spectra a round the clam i ed O VII lines at z = 0.027 and z = 0.01 1, respectively. The 
red triangle and blue diamond indicate the detections and their la ranges reported by Nicastro et al. 12005). In simulations, the instrumental response files were 
taken from Chandra observation with ID 4148 and XMM-Newton observation with ID 0099280201. 

ligible and the K/3 line becomes too weak to be detectable. 
Therefore, we have to rely on the instrumental resolution 
(R ~ 3600) to resolve the line and obtain the dynamical and 
kinematic information of O vn-absorbers. Please note that, to 
date, the X-ray spectrograph with the highest spectral resolu- 
tion is the Hig h Energy Trans mission Grating Spectrograph 
(HETGS: ICanizares et al.Nl 9861) aboard on Chandra, which 
has two sets of gratings: the high energy grating (HEG) and 
medium energy grating (MEG) with resolutions of ~ 1700 
and 900 at 20 A. Clearly, they are still not high enough for 
resolving the expected WHIM absorption line. Furthermore, 
the small effective area (< 1 cm 2 and ~ 2 cm 2 ) around 20 
A makes the HETGS infeasible for the task discussed here 
(please also see below). Figure [9] shows the necessity for 
high spectral resolution by comparing the LSFs of XMM- 
Newton, Chandra, and a next-generation X-ray spectrograph 
with R = 4000 and the corresponding simulated spectra. It 
demonstrates that an absorber with two velocity components 
separated by 150 km s" 1 cannot be resolved by XMM-Newton 
(Chandra) even from a spectrum with S/N as high as 350 
(200), but it can be clearly resolved by the next-generation 
spectrograph from a spectrum with S/N ~ 30. 

With R = 4000, an effective area of A > 100 cm 2 is needed 
for the spectrograph to systematically survey the WHIM O vn 
along 15 sight lines within 6 months. Assuming that all O vn 
absorbers have b v = 75 km s" 1 , it would take 11.2 Ms for a 
spectrograph with R = 4000 and A =100 cm 2 t o fini sh survey- 
ing those 15 QSOs with >0.2 mCrab (Section |4~TT i and mea- 
sure the intervening WHIM absorbers with Afovn ~ 10 15 cm -2 
at > 4cr significance. Assuming a 70% on-target timing ef- 
ficiency, this is equivalent to a half year. To constrain the 
absorbers with jVovii ~ 10 14 5 cm -2 with the same amount of 
exposure time, A = 1000 cm 2 is required (Figure[8]). 

In this section, we focused on the instrumental requirement 
particular for the study of the WHIM. In fact, X-ray absorp- 
tion/emission line diagnostics can also provide essential in- 
formation for the physics of the ISM, stellar coronae, black 
hole accretion and outflow, quasar accretion and feedback, 
etc., and remarkable progress has already be en made with 
Chandra and XMM-Newton observations (e.g.. [Yao & Wang I 
l2006t ITesta et alJl2008t lMiTIerl^l2T)06t fArav et al J 120071) 7 
With an X-ray spectrograph equipped with R ~ 4000 and 
A ~ 100 cm 2 , we should be able to explore all these fields 
in great detail. 



NOTE. — Required exposure times in unit of ks for future X-ray spectro- 
graphs with various spectral resolutions to obtain > 4<r significance measure- 
ment of O VII Ka with Afovil = 10 15 cm -2 and different dispersion velocities. 
The values in parentheses are required exposures for > 6a constraint (Fig- 
ure[8j. 

with 0.2 mCrab flux. 

Our simulations results deserve some explanation. First, 
each significance contour has a steep PL plus a relatively flat 
floor in the logarithm FAE-Afovii space, which roughly cor- 
responds to the shape of the COG (Figure [6j. The transition 
locations indicate the column densities from which the satu- 
ration of a line begins to affect the measurement. Second, to 
constrain Afovii = 10 15 cm" 2 at the same significance levels, it 
takes an increasingly long time with increasing b v for all res- 
olutions except for R = 1500 (Tabled. This is mainly due to 
line broadening. With R = 1500, no lines with b v < 120 km s" 1 
can be resolved, since the intrinsic line width is narrower than 
the instrumental LSF Therefore the simulations are analogous 
to those for Chandra and XMM-Newton, which explains the 
required exposure times being essentially same for all b v val- 
ues. At /?>3600, a line with b v > 50 km s" 1 becomes re- 
solved. Because the column density Novn = 10 15 cm -2 lies on 
the linear part of the COG (Figure |6j, an increased b v does 
not increase the EW but broadens the line profile, and there- 
fore requires longer exposure to pick up the absorption signals 
from the continuum. 

A spectrograph with spectral resolution of R > 3600 is re- 
quired to resolve the expec ted W HIM O vn absorption lines. 
As we discussed in Section |4~T1 the lower limit to the disper- 
sion velocity of the WHIM O vn absorbers is expected to be 
£>„>50km s , which is equivalent to aFWHM of 83.3 km s" 1 
for a Gaussian profile fit. Although we can measure b v by 
comparing the different saturation levels of the Ka and K/3 
transitions of O vn from the "low" resolutio n observations 
made with Chandra and XMM-Newton (e.g., lYao & Wang I 
120051 120071) . such technique cannot be applied to a column 
density as low as Afovii ~ 10 15 cm" 2 because saturation is neg- 
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FIG. 8. — The required exposures for detecting O VII Ka line at 2, 4, 6 and 8<r significance levels (corresponding to curves from bottom to top in each panel) 
for the absorbing gas with dispersion velocities of 50, 75, and 100 km s _1 (panels from left to right) and column densities logAfovufcm -2 ) = 14.5 — 17.0. The top 
X-axis shows the corresponding equivalent width in mA. The simulations used flux integrated from 0.5-2 keV and assumed a Gaussian profile of the line spread 
function of a spectrograph with spectral resolution R = A/FWHM = 1500, 3000, 4000, 6000 (panels from top to bottom), respectively. 



5. SUMMARY 

1. We analyzed all available (as of 2011 April) grating 
data and obtained two coadded spectra with S/N~90 per 
50-mA resolution element from Chandra observations, and 
S/N~190 from XMM-Newton observations. Neither Chan- 
dra nor XMM-Newton observations support t he existence of 
the tw o WHIM systems previously reported by Nic astro et al. 
(2005). 

2. We ran bootstrap simulations for the detecting limits of 
the current X-ray telescopes. We find that the reported EWs 
of the O vii Ka at z = 0.01 1 and z = 0.027 should have been 
measured at > 3 J a and > 2.3c, in contrast to the fact that we 
only obtained upper limits to the EWs. 

3. According to the numerical simulations, the WHIM ab- 
sorbers with A'ovii ~ 10 15 cm" 2 could sample > 30-50% of the 
O vn-bearing baryons. To systematically survey the 15 QSO 
sight lines along which the IGM O vi absorbers have been de- 
tected, future X-ray telescopes should be able to facilitate the 



WHIM study via the X-ray absorption line spectroscopy from 
background QSOs with fluxes of ~ 0.2 mCrab to find > 10 
WHIM systems. It takes impractical long (341 Ms for Chan- 
dra and 804 Ms for XMM-Newton) exposure time for current 
X-ray telescopes to complete the survey. 

4. It would require ^11 Ms on-target exposures for fu- 
ture X-ray spectrographs equipped with spectral resolution 
R > 4000 and effective area A > 100 cm 2 to finish the survey. 

We are grateful to Britton Smith for making his numeri- 
cal simulation results available to us before publication, and 
to Charles Danforth for extensive discussion on this work. 
We thank the anonymous referee for the constructive com- 
ments that help to greatly clarify the manuscript. This work 
is supported by NASA through ADP grant NNX10AE86G, 
and grant NNX08AC14G to the University of Colorado for 
data analysis and scientific discoveries related to the Cosmic 
Origins Spectrograph on the Hubble Space Telescope. 
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